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ABSTRACT:  Normal  and  surface-enhanced  Raman  spectra  for  a  set  of  substituted 
benzenethiols  were  measured  experimentally  and  calculated  from  static  polarizability 
derivatives  determined  with  time-dependent  density  functional  theory  (TDDFT).  Both  silver 
and  gold  cluster— thiolate  complexes  were  studied  to  investigate  how  the  chemical 
enhancement  varies  with  substituent.  The  experimental  relative  peak  intensities  and  positions 
are  well-matched  by  their  theoretical  counterparts.  The  static  chemical  enhancement  of  the 
ring  stretching  modes  near  1600  cm-1  is  determined  experimentally  and  computationally  for 
each  derivative,  and  it  is  found  that  the  experimental  enhancement  varies  by  a  factor  of  10  as  a 
result  of  chemical  substitution,  with  stronger  electron  donating  groups  on  the  benzene  unit 
leading  to  higher  enhancements.  The  calculated  trends  with  substitution  match  experiment 
well,  suggesting  that  TDDFT  is  describing  the  chemical  effect  qualitatively,  if  not 
quantitatively,  in  the  static  (low-frequency)  limit.  A  two-state  model  is  developed,  providing 

qualitative  insight  into  the  results  in  terms  of  the  variation  of  ligand-to-metal  charge-transfer  excitation  energy  with  substitution. 
SECTION:  Spectroscopy,  Photochemistry,  and  Excited  States 


Raman  spectroscopy  induces  and  interrogates  inelastic 
scattering  of  light  by  molecules,  with  the  scattered  light 
providing  a  vibrational  fingerprint  unique  to  each  molecule. 
This  inherently  weak  process  has  been  shown  to  be  greatly 
enhanced  when  molecules  are  adsorbed  onto  noble-metal 
surfaces.1  This  leads  to  surface-enhanced  Raman  scattering 
(SERS),  a  powerful  analytical  tool.  Signal  enhancement  is 
widely  attributed  in  literature2-8  to  two  main  mechanisms:  the 
chemical  enhancement  (CE)  and  electromagnetic  (EM) 
mechanisms. 

The  EM6,7  involves  amplification  of  the  local  electric  field  at 
the  molecule’s  position  due  to  plasmon  excitation  of  the 
metallic  structure.  Enhancements  in  the  range  of  104  to  108  are 
commonly  reported.  There  is  good  evidence  that  this 
mechanism  governs  the  dependence  of  the  SERS  signal  on 
nanoparticle  morphology.  The  CE  mechanism  is  not  as  well 
understood,  but  it  is  postulated  to  arise  from  three  distinct 
processes:8-10  enhancements  from  molecular  excitation  reso¬ 
nances,  charge-transfer  resonances,  and  nonresonant  changes  in 
the  molecular  polarizability  that  occur  upon  adsorption  onto 
the  metal  surface.  There  have  been  attempts  to  develop 
theories  that  include  all  contributions  of  the  SERS  effect 
together,  such  as  work  by  Lombardi  et  al.  based  on  Herzberg— 
Teller  relations.9,11  However,  investigating  these  effects 
separately  allows  for  the  magnitude  of  each  contribution  to 
be  isolated  and  quantified. 

Herein  we  focus  on  the  third  CE  mechanism,  which  is  also 
called  the  static  chemical  enhancement  (CHEM).  This 
enhancement  mechanism  is  intrinsically  linked  with  the 
creation  of  a  new  chemical  system  (i.e.,  the  metal— molecule 


complex)  leading  to  changes  in  the  polarizability  and 
vibrational  modes  of  the  molecular  moeity.  Changes  in  the 
molecular  polarizability  upon  binding  are  believed  to  lead  to 
enhancements  in  the  range  of  ~10°  to  103  for  Raman 
scattering,5’12'13  but  this  estimate  carries  uncertainty  due  to 
the  lack  of  careful  theory/ experiment  comparisons.  Recently, 
Morton  and  Jensen  have  proposed  a  two-state  model  for 
studying  the  static  CHEM  enhancement.14  They  used  TDDFT 
calculations  in  the  static  (low-frequency)  limit  for  molecules 
adsorbed  onto  silver  clusters  to  define  the  chemical  enhance¬ 
ment  factor  (through  comparison  with  the  calculated  intensity 
in  the  absence  of  surfaces).  They  found  that  molecule— metal 
systems  that  stabilize  pi  backbonding  provide  an  enhancement 
to  the  attached  molecule’s  polarizability  and  hence  to  enhanced 
Raman  scattering.  They  argued  that  calculations  of  this  effect 
for  small  silver  clusters  would  give  results  similar  to  realistic 
SERS  nanoparticle  substrates,  as  the  static  limit  Raman 
intensity  is  only  weakly  dependent  on  cluster  size.  The 
chemical  enhancement  factor  they  calculated  was  purely 
theoretical;  however,  comparisons  of  the  TDDFT  results  with 
a  two-state  model  were  provided  to  reveal  the  analytical  origin 
of  the  effect.  This  comparison  showed  that  the  enhancement 
factor  typically  depends  on  the  energy  gap  between  the  HOMO 
of  the  cluster  and  the  LUMO  of  the  adsorbed  molecule,  as  is 
physically  plausible.  However  a  concern  with  this  formulation  is 


Received:  June  14,  2013 
Accepted:  July  22,  2013 


ACS  Publications  ®  xxxx  American  Chemical  Society 


2599 


dx.doi.org/10.1021/jz4012383  \J.  Phys.  Chem.  Lett.  2013,  4,  2599-2604 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

2013 


2.  REPORT  TYPE 


4.  TITLE  AND  SUBTITLE 

A  Look  at  the  Origin  and  Magnitude  of  the  Chemical  Contribution  to  the 
Enhancement  Mechanism  of  Surface-Enhanced  Raman  Spectroscopy 
(SERS):  Theory  and  Experiment 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Northwestern  University, Department  of 
Chemistry  ,Evanston,IL, 60208-3113 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED 

00-00-2013  to  00-00-2013 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

Normal  and  surface-enhanced  Raman  spectra  for  a  set  of  substituted  benzenethiols  were  measured 
experimentally  and  calculated  from  static  polarizability  derivatives  determined  with  time-dependent 
density  functional  theory  (TDDFT).  Both  silver  and  gold  cluster &#8722;thiolate  complexes  were  studied  to 
investigate  how  the  chemical  enhancement  varies  with  substituent.  The  experimental  relative  peak 
intensities  and  positions  are  well-matched  by  their  theoretical  counterparts.  The  static  chemical 
enhancement  of  the  ring  stretching  modes  near  1600  cm&#8722;l  is  determined  experimentally  and 
computationally  for  each  derivative,  and  it  is  found  that  the  experimental  enhancement  varies  by  a  factor 
of  10  as  a  result  of  chemical  substitution,  with  stronger  electron  donating  groups  on  the  benzene  unit 
leading  to  higher  enhancements.  The  calculated  trends  with  substitution  match  experiment  well,  suggesting 
that  TDDFT  is  describing  the  chemical  elect  qualitatively,  if  not  quantitatively,  in  the  static 
(low-frequency)  limit.  A  two-state  model  is  developed,  providing  qualitative  insight  into  the  results  in  terms 
of  the  variation  of  ligand-to-metal  charge-transfer  excitation  energy  with  substitution. 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

6 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Letter 


The  Journal  of  Physical  Chemistry  Letters 


Figure  1.  Left  pane:  Structure  of  the  Au19— benzenethiolate  complex  used  in  TDDFT  calculations.  The  metal  cluster  is  constructed  by  truncating  one 
vertex  of  a  20  atom  tetrahedral  cluster.  Right  pane:  Surface  of  gold  INRA  SERS  substrate  used  for  Raman  experiments.  The  substrates  were 
fabricated  by  depositing  200  nm  of  metal  over  600  nm  silica  microspheres. 


Table  1.  Vibrational  Mode  Energies  for  Benzenethiol  and  Benzenethiol— Metal  Complexes" 


BT  (cm”1) 

Ag-BT  (cm  1 ) 

Au-BT  (cm  *) 

exp.  BT  (cm  1 ) 

exp.  Ag-BT  (cm  1 ) 

exp.  Au-BT  (cm  1 ) 

mode  description 

689.6 

682.9 

683.4 

699 

688 

688 

ring  and  C— S  stretch 

989.3 

987.6 

988.3 

1000 

996 

996 

ring  breathing 

1017.5 

1013.1 

1013.3 

1025 

1018 

1018 

ring  breathing 

1073.6 

1079.9 

1054.3 

1054.1 

1092 

1068 

1069 

ring  breathing  and  C— S  stretch 

1578.0 

1564.7 

1565.1 

1584 

1573 

1573 

C— C  stretch 

“First  three  columns  on  the  left  are  from  TDDFT  calculations,  while  columns  denoted  “exp”  represent  modes  observed  experimentally. 
Experimental  SERS  spectra  were  recorded  using  metal  INRA  SERS  substrates. 


that  pure  density  functionals  tend  to  underestimate  energy  gaps 
between  occupied  and  virtual  states,  so  it  was  not  clear  that  the 
chemical  effects  were  correctly  described.  We  show,  with 
extensive  comparisons  to  experiment,  that  the  same  level  of 
theory  (TDDFT  in  the  static  limit)  can  reliably  calculate  the 
contribution  of  adsorption-induced  nonresonant  molecular 
polarizability  changes  to  the  SERS  effect. 

Substrates  used  in  SERS  studies  usually  consist  of  Ag  or  Au 
nanoparticles7'15  or  roughened  surfaces.16  18  In  general,  such 
substrates  are  designed  to  provide  large  EM  enhancement, 
albeit  with  considerable  spot-to-spot  signal  variation  across  the 
surface.  For  this  work,  it  is  important  to  have  uniform  EM 
enhancement  over  large  areas,19’20  as  this  simplifies  the 
separation  of  EM  and  CE.  Both  silver  and  gold  immobilized 
nanorod  assembly  (INRA)21,22  substrates  have  been  shown  to 
possess  these  traits,  leading  to  their  use  in  this  work. 

Theoretical  static  SERS  spectra  including  chemical  effects 
have  been  generated  by  incorporating  a  small  metal  cluster  into 
the  TDDFT  calculations.  This  follows  our  previous  work, 
where  a  20  atom  cluster  with  approximately  tetrahedral  shape 
was  found  to  provide  a  reasonable  model  of  a  (ill)  surface 
(any  of  the  faces  of  the  tetrahedron)  upon  which  a  molecule  is 
adsorbed.1”'13  In  the  present  application,  this  choice  is 
complicated  by  the  thiol  molecule’s  propensity  to  undergo 
chemisorption  on  these  surfaces.  After  studying  several 
alternatives,  we  chose  a  Ag19/Au19— thiolate  complex,  with  the 
vertex  atom  opposite  the  adsorption  site  removed,  as  this 
closed-shell  structure  provides  spectra  with  the  best  exper¬ 
imental  agreement.  A  representative  structure  of  modeled 
molecule— metal  complexes  is  shown  in  Figure  1.  Geometry 
optimizations  were  performed  utilizing  a  starting  structure  with 
the  sulfur  above  the  central  metal  atom  on  the  triangular  cluster 
face.  Local  minima  correspond  to  structures  where  sulfur 
bridges  two  metal  atoms.  For  the  gold  clusters,  some  of  the 
benzenethiol  derivatives  also  have  a  local  minimum  structure 
with  sulfur  bound  to  the  central  metal  atom  of  the  face. 


However,  these  structures  were  always  higher  in  energy  than 
the  bridging  structures.  The  bridging  structures  were  therefore 
chosen  as  the  metal— molecule  complex  used  to  calculate  SERS 
spectra.  Note:  this  bridging  structure  is  related  to  S— Au— S 
bridges  thought  to  be  important  for  arylthiolates  adsorbed  onto 
small  gold  clusters23  26  and  on  Au(lll)  surfaces.27 

Seven  molecules  were  studied  experimentally  and  theoret¬ 
ically  using  both  gold  and  silver  substrates.  Not  all  molecules 
were  stable  on  the  metal  substrates,  leading  to  a  smaller 
experimental  data  set  with  eight  molecule/surface  combina¬ 
tions.  Chemical  enhancement  factors  (EFs)  were  determined 
for  three  molecules  on  silver  surfaces:  benzenethiol  (BT),  4- 
mercaptophenol  (MP),  and  4-mercaptobenzoic  acid  (MBA) 
and  five  on  gold  surfaces:  BT,  MBA,  4-methoxybenzenethiol 
(MBT),  4-ethoxybezenethiol  (EBT),  and  4-nitrobenzenethiol 
(NBT). 

Calculated  energies  of  the  benzenethiol  derivatives’  vibra¬ 
tional  modes  are  shifted  when  the  molecules  are  bound  to  the 
clusters.  Select  normal  mode  energies  are  presented  in  Table  1 
for  benzenethiol  in  the  gas  phase  and  bound  to  silver  and  gold 
clusters.  Large  shifts  occur  for  the  ring-breathing  mode  near 
1060  cm-1  and  the  ring  C— C  stretches  near  1570  cm-1,  while 
the  modes  at  ~6 12  and  988  cm-1  see  minor  changes.  Cartesian 
coordinates  and  vibrational  normal  modes  between  300  and 
2000  cm-  for  all  optimized  structures  and  detailed 
comparisons  of  calculated  and  experimental  spectra  can  be 
found  in  the  Supporting  Information.  To  further  verify  our 
theoretical  description  of  normal  and  SERS  spectra  as  well  as 
the  static  chemical  effects  involved,  CHEM  EFs  for  several 
vibrational  modes  of  each  molecule  were  determined.  A  linear 
trend  is  expected  when  comparing  theoretical  and  experimental 
results  if  the  theoretical  model  is  a  reasonable  representation  of 
the  experimental  system.  This  technique  has  been  previously 
used  by  Zayak  et  al.  to  determine  surface  binding  geometries.28 
This  analysis  is  also  independent  of  the  surface  coverage  of  the 
molecules  (as  long  as  there  is  no  change  in  binding  geometry). 
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Peaks  were  chosen  for  analysis  if  they  were  present  in  both  the 
normal  and  SERS  experimental  spectra.  As  expected,  some 
peaks  in  the  solution  Raman  spectra  were  obscured  by  solvent. 

Vibrational  modes  that  could  be  properly  analyzed  are 
plotted  in  Figure  2.  For  clarity,  EFs  for  each  peak  are  scaled 


Figure  2.  Comparison  of  scaled  experimental  and  simulated 
enhancement  factors  for  various  modes  of  the  molecule  and  choice 
of  substrate.  Circles  denote  Ag  substrate,  squares  denote  Au  substrate. 
Colors  denote  the  substituted  benzenethiol:  BT  (black),  NBT  (red), 
MBT  (green),  EBT  (light  blue),  MP  (blue),  MA  (purple).  The 
straight  line  is  y  =  x.  Experimental  SERS  spectra  were  recorded  using 
metal  INRA  SERS  substrates. 


such  that  the  largest  EF  for  each  molecule  is  set  in  the  range  0.8 
to  1.0.  Values  for  NBT  and  MA  on  silver  are  not  included 
because  these  molecules  undergo  an  induced  photoreac¬ 
tion.29'30  Only  a  single  peak  (~1600  cm-1)  was  sufficiently 
differentiated  from  solvent  for  MBA  on  both  metals.  Of  the 
remaining  combinations,  the  MP  on  gold  and  the  MBT  and 
EBT  on  silver  show  very  poor  agreement  of  calculated  and 
experimental  chemical  EFs.  For  these  complexes,  it  is  possible 
that  molecules  were  bound  to  the  surface  in  a  geometry 
different  from  what  was  assumed.  Also,  extensive  delamination 
was  observed  for  EBT  on  the  silver  INRA  substrates.  It  is  likely 
that  the  MP  and  MBT  are  delaminating  similarly  on  a  smaller 
scale.  Because  this  behavior  results  in  unreliable  EFs,  these 
combinations  were  not  considered  in  further  analysis. 
Comparisons  involving  the  remaining  molecules  result  in  linear 
behavior,  with  slopes  near  unity  and  intercepts  near  zero.  The 
EF  calculation31  used  for  this  study  is  described  in  detail  in  the 
Supporting  Information. 

Chemical  EFs  for  the  set  of  peaks  occurring  at  ~1600  cm-1 
for  each  of  the  benzenethiol  derivatives  determined  from 
experiment  and  TDDFT  calculations  are  shown  in  Table  2. 
Peaks  at  ~  1 600  cm-1  were  chosen  because  they  are  clear  of 
interference  (including  solvent  peaks)  in  the  experimental 
SERS  and  normal  Raman  spectra  for  the  entire  molecule  set. 
The  spectra  corresponding  to  the  silver  SERS  measurements 
are  included  in  Figure  3.  Vibrational  modes  of  the  two  peaks 
near  1600  cm-1  included  in  the  EF  calculations  are  the  two 
nearly  degenerate  ring-stretching  modes  shown  in  Figure  4. 
The  position  and  splitting  of  these  modes  are  dependent  on  the 
identity  of  the  substituent.  Whether  one  or  two  peaks  were 
observed  experimentally  in  this  spectral  region,  the  entire  area 
under  the  curve  was  used  to  calculate  EFs. 


Table  2.  Chemical  Enhancement  Factors  Determined  from 
Experiments  and  Calculated  Using  TDDFT" 


experimental  1600  cm  1 
EF  (scaled) 

TDDFT  1600  cnT1  EF 
(absolute) 

Au- 

-oc2hs 

53.1 

111.0 

-och3 

45.5 

109.8 

-H 

36.3 

36.3 

-co2h 

14.7 

15.7 

-no2 

14.9 

19.4 

Ag- 

-OH 

60.5 

60.6 

-H 

11.4 

11.4 

-co2h 

5.9 

7.8 

"Experimental  EF  values  for  each  metal  were  scaled  relative  to  the 
benzenethiol  result,  with  BT  being  scaled  equal  to  the  corresponding 
theoretical  value.  Experimental  SERS  spectra  were  recorded  using 
metal  INRA  SERS  substrates. 


350  600  850  1100  1350  1600 

Raman  Shift  (cm-1) 


Figure  3.  SERS  spectra  of  benzenethiol  derivatives  on  AgINRA 
surface.  All  measurements  were  made  with  a  785  nm  excitation  source. 
The  EF  1600  values  reported  correspond  to  the  enhancements  calculated 
using  the  singlet  or  doublet  of  peaks  found  ~1600  cm-1. 


Figure  4.  Vibrational  motions  of  benzenethiol  near  1600  cm  1  chosen 
for  EF  comparison  in  Table  2  and  Figure  5. 


Effects  due  to  molecular  resonances  were  avoided  by  using 
incident  light  (785  nm)  not  energetic  enough  to  electronically 
excite  any  states  associated  with  the  molecules  considered.32 
The  enhancement  observed  is  therefore  a  combination  of  the 
EM  enhancement  from  the  plasmonic  substrate  and  a  static 
chemical  effect  due  to  ground-state  molecule— surface  inter¬ 
actions.  Dividing  the  SERS  intensity  per  molecule  by  the 
normal  Raman  intensity  per  molecule  gives  an  EF  measurement 
that  is  composed  only  of  the  EM  enhancement  (from  the 
substrate)  and  of  static  chemical  effects.  The  uniformity21  of  the 
INRA  substrate  assures  that  the  EM  enhancement  is  a  constant 
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contributor.  Thus  the  variation  in  measured  EF  for  different 
molecules  is  a  direct  reflection  of  the  static  chemical 
enhancement.  Theoretical  calculations  of  the  Raman  spectra 
were  performed  in  the  static  limit  to  determine  Raman  cross 
sections  for  both  the  free  thiol  and  the  corresponding  thiolate— 
metal  cluster  complex. 

Figure  5  presents  the  comparison  of  experimental  and 
simulated  EFs  for  the  1600  cm-1  peaks.  In  contrast  with  Figure 


0.0  0.2  0.4  0.6  0.8  1.0 

Scaled  Theoretical  EF 

Figure  5.  Comparison  of  scaled  experimental  and  simulated 
enhancement  factors  for  the  1600  cm-1  peak(s).  Ag  complexes  are 
red  points.  Au  complexes  are  blue  points.  The  gray  line  is  y  -  x. 
Experimental  SERS  spectra  were  recorded  using  metal  INRA  SERS 
substrates. 


2,  this  plot  shows  a  more  linear  relationship  between 

experiment  and  theory,  with  both  data  showing  an  increase  in 

enhancement  as  the  substituent  group  becomes  more  electron- 

donating.  Linear  fitting  of  the  data  gives  a  slope  near  unity, 

suggesting  that  TDDFT  is  reasonably  accurate  in  calculation  of 
relative  chemical  EFs. 

The  absolute  CHEM  EFs  are  calculated  to  be  between  101 

and  102 * * *,  which  is  in  the  range  expected  based  on  previous 
work.14  On  gold,  the  relative  enhancements  from  highest  to 
lowest  vary  by  a  factor  of  3.6  in  the  experimental  measurements 
and  a  factor  of  7  in  the  density  functional  calculations.  On 
silver,  the  factors  are  10.3  and  8,  respectively.  In  agreement  with 
other  studies,33  the  absolute  CHEM  enhancements  were 
calculated  to  be  larger  on  gold  than  on  silver.  This  is  due  to 
the  lower  energy  of  the  LUMO  of  the  gold  cluster  than  that  of 


the  silver  cluster,  as  previously  noted  by  Morton  et  al.8  Both 
theory  and  experiment  agree  that  the  electron-donating 
substituents  give  rise  to  the  largest  CHEM  enhancement. 

Now  let  us  consider  comparisons  of  the  experimental  and 
theory  results  with  those  from  the  two-state  model  developed 
by  Morton  and  Jensen.  Table  3  displays  the  comparison  of 
experimental  and  TDDFT  integrated  EFs  with  two-state  model 
results.  Assuming  as  was  done  in  previous  studies  of  pyridine 
derivatives  that  the  important  transition  is  the  metal  HOMO  to 
molecule  LUMO  (metal-to-ligand  charge  transfer  or  “MLCT”), 
the  two-state  model  predicts  a  reduction  in  intensity  that  is 
unphysical  in  nature.  Considering  instead  a  molecular  HOMO 
to  metal  LUMO  transition8  (ligand-to-metal  charge  transfer  or 
“LMCT”),  the  experimental,  TDDFT,  and  two-state  model 
results  all  display  the  same  trend.  Note:  Exact  agreement  is  not 
expected  with  the  approximations  made  in  this  simple  model. 
In  addition,  there  are  uncertainties  in  how  the  experimental  and 
TDDFT  EFs  are  determined  depending  on  which  vibrational 
modes  are  chosen  for  comparison. 

For  the  experimental  EFs  only,  the  observed  peaks  for  each 
molecule  rather  than  all  vibrational  modes  were  selected  for 
comparison.  In  the  TDDFT  results,  we  can  either  choose  the 
same  modes  as  in  the  experiment  or  we  can  integrate  over  all 
modes.  To  see  the  size  of  the  difference  between  the  results, 
TDDFT-integrated  EFs  were  calculated  using  both  the  full 
range  of  peaks  between  500  and  2000  cm-1  and  the  subset 
(labeled  "partial”  in  Table  3)  used  for  the  experimental 
determination.  The  results  show  variations  by  a  factor  of  2  for 
some  substituents.  The  comparisons  between  experiment  and 
TDDFT  are  generally  improved  by  the  partial  selection 
procedure.  The  two-state  model  shows  trends  with  substituents 
that  are  greatly  exaggerated  compared  with  either  experiment 
or  TDDFT  but  are  qualitatively  similar.  We  must  emphasize 
that  the  two-state  model  is  a  simplification  and  in  the  present 
application  we  have  used  a  and  b  factors  from  the  study  by 
Morton  and  Jensen,  which  assumed  a  MLCT  mechanism.  It  is 
not  clear  from  our  work  if  revising  these  parameters  for  the 
LMCT  version  of  the  model  would  improve  comparisons  with 
TDDFT. 

To  summarize,  it  was  shown  through  comparison  with 
experiment  that  TDDFT  with  a  BT-Ag19  or  BT-Au19  “SERS” 
model  accurately  predicts  chemical  EFs  across  a  range  of 
substituted  benzenethiols.  The  calculations  provided  a 
satisfactory  prediction  of  relative  chemical  enhancements  for 
peaks  within  the  same  spectrum.  These  calculations  show  that 
there  is  a  modest  CHEM  enhancement  in  the  range  of  ~10°  to 


Table  3.  Integrated  Enhancement  Factors  Determined  from  Experiments,  Calculated  Using  TDDFT,  and  Calculated  Using  a 
Two-State  Model  Considering  Two  Different  Sets  of  States" 


experimental  integrated  EF  (x  106 *) 

TDDFT  integrated  EF 

TDDFT  partial  integrated  EF 

MLCT  two-state  EF 

LMCT  two-state  EF 

-oc2H5 

66 

51 

140 

0.03 

780 

-och3 

38 

54 

66 

0.03 

900 

-H 

31 

24 

28 

0.12 

88 

-co2h 

10 

15 

16 

0.13 

7.4 

-no2 

23 

13 

25 

0.18 

1.3 

-OH 

58 

27 

46 

0.02 

24 

-H 

20 

10 

12 

0.04 

29 

-co2h 

19 

8.0 

7.8 

0.54 

0.25 

“Experimental  SERS  spectra  were  recorded  using  metal  INRA  SERS  substrates. 
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102  for  substituted  benzenethiols,  with  more  strongly  donating 
groups  leading  to  larger  enhancements.  The  two-state  model 
proved  to  be  useful  for  “rough”  comparison  between  theory  and 
experiment;  however,  we  found  that  it  was  necessary  to  use  the 
LMCT  energy  gap  to  determine  meaningful  EFs  from  this 
model,  while  using  the  MLCT  gap  results  even  in  qualitative 
errors. 

■  EXPERIMENTAL  METHODS 

BT,  MP,  MBT,  NBT,  and  MBA  were  used  as  received  from 
Sigma-Aldrich.  EBT  was  purchased  from  Santa  Cruz  Bio¬ 
technology.  All  solvents  were  purchased  from  Sigma-Aldrich. 
BT,  MP,  MBT,  and  EBT  were  used  for  overnight  incubation  at 
150  mM  in  ethanol.  NBT  and  MBA  were  prepared  at  the  same 
concentration  in  acetone  and  isopropanol,  respectively. 

Metal  INRA  SERS  substrates  were  fabricated  as  previously 
described.21'22  Gold  and  silver  films  (200  nm  thick)  were 
deposited  at  a  rate  of  1  A-s_1  under  high  vacuum  (2.0  X  10-7 
Torr)  over  the  nanosphere-covered  surface  using  thermal  vapor 
deposition  (Kurt  J.  Lesker,  PVD  75).  The  metal  mass  thickness 
and  deposition  rate  were  measured  by  a  6  MHz  gold-plated 
quartz  crystal  microbalance  (Sigma  Instruments,  SQM-160). 
The  gold  INRA  surfaces  were  plasma-cleaned  for  5  min  with 
99.999%  02  in  a  refractive  ion  etcher  (South  Bay  Technology, 
RIE-2000)  at  196  mT  and  51  W.  An  inverted  microscope 
(Nikon  Eclipse  Ti)  was  used  for  all  Raman  measurements.  An 
Innovative  Photonics  Solutions  stabilized  785  nm  diode  laser 
was  used  for  single-frequency  experiments.  The  detection 
system  was  an  imaging  spectrograph  (Acton  SpectraPro  2300i) 
with  an  LN2-cooled  back-thinned  deep-depletion  CCD 
detector  (Roper,  model  spec  10:400  BR,  1340  X  400  pixels). 
The  source  was  directed  through  the  back  entrance  of  the 
microscope  and  focused  down  to  a  70  //m  X  70  pm  spot  onto 
the  opaque  metal  INRA  substrates  using  a  20X  microscope 
objective  (wd  =  2.1  mm). 

Computational  Methods.  Calculations  presented  in  this  work 
were  performed  using  the  Amsterdam  Density  Functional 
(ADF)  program  package.34  Full  geometry  optimization  and 
frequency  calculations  for  isolated  benzenethiols  and  the 
molecule— metal  cluster  complexes  were  completed.  A 
gradient-corrected  Becke— Perdew  (BP86)  exchange— correla¬ 
tion  functional  was  used  with  a  triple-^  polarized  Slater-type 
(TZP)  basis  set  for  all  atoms.  Polarizabilities  were  calculated 
using  the  AOResponse  module  with  the  same  level  of  theory. 
Scalar  relativistic  effects  were  accounted  for  with  the  zeroth 
order  regular  approximation  (ZORA).  Raman  intensities  were 
calculated  from  static  polarizabilities,  as  outlined  in  the 
Supporting  Information. 

■  ASSOCIATED  CONTENT 
G  Supporting  Information 

Optimized  geometries,  vibrational  frequencies,  and  theoretical 
and  experimental  Raman  spectra  of  all  molecules.  Example 
experimental  Raman  and  SERS  measurements  along  with  an  EF 
calculation  and  further  theory  computational  detail.  This 
material  is  available  free  of  charge  via  the  Internet  at  http:// 
pubs.acs.org. 
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